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ABSTRACT 



We compute, for several solar models, the rates P at 
which the solar radial p modes are expected to be ex- 
cited. The solar models are computed with two differ- 
ent local treatments of convection : the classical mixing- 
length theory (MLT hereafter) and ICanuto et all jl99(Tl 
CGM hereafter)'s formulation. 

For one set of solar models (EMLT and ECGM models), 
the atmosphere is gray and assumes Eddington's approx- 
imation. The models only assume one mixing-length pa- 
rameter and reproduce the solar radius at solar age but 
not the Balmer line profiles. For a second set of mod- 
els (KMLT and KCGM models), the atmosphere is built 
using a T(r) law which has been obtained from a Ku- 
rucz's model atmosphere computed with the same local 
treatment of convection. The mixing-length parameter in 
the model atmosphere is chosen so as to provide a good 
agreement between synthetic and observed Balmer line 
profiles, while the mixing-length parameter in the inte- 
rior model is calibrated so that the model reproduces the 
solar radius at solar age. 

For the MLT treatment, the rates P do depend signifi- 
cantly on the properties of the atmosphere. Indeed dif- 
ferences in P between the EMLT model and the KMLT 
are found very large. On the other hand, for the CGM 
treatment, differences in P between the ECGM and the 
KCGM models are very small compared to the error bars 
attached to the seismic measurements. The excitation 
rates P for modes from the EMLT model are significantly 
under-estimated compared with the solar seismic con- 
straints. The KMLT model results in intermediate values 
for P and shows also an important discontinuity in the 
temperature gradient and the convective velocity. On the 
other hand, the KCGM model and the ECGM model yield 
values for P closer to the seismic data than the EMLT and 
KMLT models. We conclude that the solar p-mode exci- 
tation rates provide valuable constraints and according to 
the present investigation cleary favor the CGM treatment 
with respect to the MLT. 

Key words: turbulence, convection, oscillations, excita- 
tion. Sun. 



1. INTRODUCTION 



We study the implication of two different local treatments 
of convection on the calculation of the rates P at which 
solar p modes are excited by turbulent convection. For 
this purpose, solar models are built with two different 
local treatments of convection : the classical MLT and 
ICanutoetanjl996 ')'s formulation (Sect.|3}. 

For a first set of models (EMLT and ECGM models), the 
atmosphere is gray and assumes Eddington's approxima- 
tion. Each of these models therefore uses one single value 
of the mixing-length parameter and reproduces the solar 
radius at solar age but not the Balmer lines. 

For a second set of solar models (KMLT and KCGM 
models), the atmosphere is built using a T(t) law which 
has been obtained from Kurucz's model atmospheres 
computed with the same local treatment of convection as 
used in the interior. The mixing-length parameter in the 
model atmosphere (aa) is chosen so as to provide a good 
agreement between synthetic and observed Balmer line 
profiles, while the mixing-length parameter in the inte- 
rior model {oLi) is calibrated so that the model reproduces 
the solar radius at solar age. 

Calculation of the excitation rates requires the computa- 
tion of the convective velocity and the convective flux. 
This is done in Sect. 13.11 and 13.21 by paying special at- 
tention to the difficulty of matching interior and atmo- 
sphere models. We finally present the computed excita- 
tion rates P for each model and compare with the solar 
data in Sect.|4l 



CALCULATION OF THE RATE AT WHICH 
SOLAR P-MODES ARE EXCITED 



The model of stochastic excitation (MSE ) used here is ba- 
sically that of ISamadi & Goupill J2001L see also Samadi 
et al. 2004, present conference). This model provides 
an expression for P(y), the rate at which a given radial 
mode with frequency vq is excited and requires a proper 
description of the structure of the very outer layers. 
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In practice the calculation of P needs the following quan- 
tities as input: 1) The mean density, the convective ve- 
locity {v) and the entropy fluctuations (s) of the solar 
model, s is directly proportional to the ratio (Fc/v)'^ 
where Fc is the convective flux. Fc and v are calculated 
as explained in Sect. |2] 2) The eigenfunctions (^) and 
their frequencies (I'o). They are computed with the adi- 
abatic code FILOU of Tran Minh & Leon ( 199^. E{k) 
the wavenumber (k) dependency of the turbulent kinetic 
energy spectrum, E{k, v). 3) The values and depth de- 
pendence of the wavenumbers fco^ the wavenumber at 
which E{k) is maximum. 4) Xk{i^), the frequency (i/) 
component of E{k, v) 

According to the results in ISamadi et alJ ( l2n03alW) ob- 
tained on the base of a solar 3D simulation the k- 
dependency of E{k, v) is approximatively reproduced 
by an analytical spectrum called 'Extended Kolmogorov 
Spectrum' (EKS) and defined in Mus ielak et alJ lll994h . 
The zy-dependency of Xk is found to be better mod- 
elled with a Lorentzian function rather than by a Gaus- 
sian function which is usually assumed for Xfc (see 
[Samadi et al. 2004, present conference). At the top of the 
superadiabatic region, it was found that ko ~ 3.6 Mm^^ 
and decreases slowly with depth inward. In this work we 
assume a constant /cq — 3.6 Mm~^ which is a good ap- 
proximation. 



3 . THE SOLAR MODELS 



The m odels are computed with the CESAM code ( iMorell 
[1997) including th e following input physicses: CEFF 
equation of state ("Christensen-Dalsgaar d et aP Il988h . 
OPAL opacities (lelesias & Rogers 1996) data comple- 
mented by Alexander & Ferguson ( 1994) data for T < 
lOfK, both set of data being given for Noels & Gre vessd 
( IT993I) solar mixture, Caughlan & Fowler ( 198^ ther- 
monuclear reaction rates, microscopic diffusion accord- 
ing to the simplified formalism of Michaud & Proffitt 
('1993') and finally Bohm-Vitense formulation of the MLT 
(jJohm - Vitense 1958) or CGM's formulation (the same 
formulation in the interior as in the atmosphere). 

ECGM and EMLT models: The atmosphere of the ECGM 
and EMLT is gray and assumes Eddington's approxima- 
tion. The mixing-length parameter a of those models (the 
same a in the interior as in the atmosphere) is adjusted in 
order to reproduce the solar luminosity and radius at the 
solar age. These models do not reproduce the Balmer line 
profile. TableQlgives the calibrated values of the mixing- 
length parameters. 



KCGM and KMLT models: The atmosphere of the 
KCGM and KMLT are restored from a T(t) -law derived 
from Kurucz's complete ATLAS9 mo del atmospheres ac - 
cording to the procedure described in Morel et al. ('1994'). 
The fit between the interior (where diffusion approxima- 
tion is valid) and atmosphere is performed in the region 
where ti < t < t2. In the interior region where r > T2, 
the temperature gradient Vi is obtained from the MLT 
or CGM formulation. In the atmospheric region where 
r < Ti, the temperature gradient Va is computed using 
the T — T law of the model atmosphere built with the 
same formulation of convection as in the interior. In the 
transition region where ri < t < T2, in order to ensure 
the continuity of the temperature gradient, V is obtained 
by a linear interpolation of Vi and Va as a function of the 
optical depth r, as follows: 

V = /3(r)Va + (l-/3(T))Vi (1) 

where /3(r) = (t2 ~ t)/{t2 — ti). In practice we found 
that ri=5 is the minimal acceptable value for ti and 
T2=20 minimizes the discontinuity between the interior 
and the atmosphere. 

Calibration of the KCGM and KMLT models: The 
mixing-length parameter of the model atmosphere 
is adjusted so as to provide the best agreement be- 
tw een synthetic and observed Balmer l ine profiles (as 
in Ivan't Veer- Menneret & Megessier 1996. for the MLT 
treatment). For all the treatments of convection the value 
Q!a — 0.5 provides the best agreement between synthetic 
and observed Balmer line profiles. The mixing-length pa- 
rameter a,; of the internal structure is adjusted in order 
that the global model reproduces simultaneously the so- 
lar radius and luminosity at the solar age. Values of the 
calibrated mixing-length parameters are given in Table|2] 



3.1. Calculation of the convective flux 

One contribution to the driving of the oscillation modes 
comes from the advection of turbulent entropy fluctua- 
tions by the turbulent movements (the so-called entropy 
source term). It scales as the square of the convective flux 

Fc can be viewed as function of a and V: Fc = /i(V, a), 
h is given by the adopted formulation of convection (here 
MLT or CGM). In the outer region (r < ti) and in the 
interior region (r > T2) Fc is computed directly from 
the function h. In the transition region, as in Eq. (1), the 
convective flux Fc can be related to the convective flux of 



Table 1. Values of the the mixing-length parameter a of 
the ECGM and EMLT models. 



Table 2. Values of the the mixing-length parameters of 
the KCGM and KMLT models: a-i (for the interior) and 
tta f/o'" the model atmosphere). 



model 



EMLT 
ECGM 



1.763 
0.91 



model 



KMLT 
KCGM 



2.42 
0.78 



0.50 
0.50 



3 




50 100 150 200 250 300 
Ro-r [km] 



Figure 1. The computed convective is plotted versus 
the optical depth t for the KMLT model ( dashed line ), the 
EMLT model (dot dashed line), the KCGM model (solid 
line) and the ECGM model (dot dot dashed line). 



the interior (F^) and the convective flux of the atmosphere 
(F^) as foUows: 

F, = \{t)F^ + {1-X{t))fI (2) 

In Eq. (2), = /i(V.,QO, F^ = /i(Va,aa), A(t) is, 
as P{t) (Eq. 1), a function to ensure the continuity of the 
convective flux. The choice for A(t), as for /3(t), is rather 
arbitrary. We assume A(t) = /3(r) for sake of simpHcity. 
Results of the calculation of Fc{t) are shown in Fig. (1). 



3.2. Calculation of the convective velocity 

The second contribution to the driving of the oscillations 
comes from the Reynolds term and is proportionnal to 
where v is the convective velocity. 

The convective velocity (v) can be viewed as a function 
of V and a: v = /(V, a), / is given by the adopted 
formulation of convection (here MLT or CGM). 

For each model, in the outer region (i.e. r < n) and in 
the inner region (r > T2), given V and a, v is computed 
as w = ./(V, a). In the transition region (ti < r < T2) 
of the KMLT and KCGM models we must deal with two 
different values of a (namely ao, a,) but v is not a linear 
function of a. We thus face the difficulty to properly de- 
fine a convective velocity consistent with F^ in this region 
(Eq. 2). 

We then proceed as follows: V and Fc are defined by 
Eq. (1) and Eq. (2) respectively. Then at fixed V and 
r, we define an equivalent mixing-length parameter, a*, 
such that Fc = ft(V,a*). Such mixing-length param- 
eter is hence variable with depth. We next compute 
V = /(V(r), q:*(t)) which then is consistent with Fc 
ofEq. (2). 
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Figure 2. The computed convective velocity v is plotted 
versus Rq ~ r where Rq is the solar radius and r the 
radius of a given layer. The legend has the same meaning 
as in Fig. (1). 

3.3. Comments 

The CGM model results in a much smaller value for 
ai than the MLT model. Furthermore for the CGM 
model, ai is found closer to aa = 0.5. In contrast, for 
the MLT treatment, the value of a, is much larger than 
aa = 0.5. These features can be explained by the fact 
that in nearly adiabatic regions convection is more ef- 
ficient in the CGM's formulation (see details in Heiter 
et al, 2002). Differences in ai between the KGCM and 
the ECGM models are relatively smaUer than between the 
KMLT and the EMLT models. 

As shown in Fig. (1) and (2), the EMLT and KMLT mod- 
els have very different convective structures. This is due 
to large differencies between ai and aa. KCGM and 
ECGM models have very similar convective structures. 
This is a consequence of the fact that the KCGM model 
results in a value of a, close to that required for the at- 
mosphere (aa = 0.5). 

For the KMLT model, there is an important discontinu- 
ity at the bottom boundary of the transition region (i.e. 
at T = 20 or Ra — r = 100 km), especially for v. In 
contrast for the KCGM the discontinuity is much less im- 
portant. These features are again directly connected with 
the differences between between ai and 

Conclusion: both the KMLT and the KCGM models re- 
produce the B aimer line profile and the solar radius and 
luminosity but the CGM model models the transition be- 
tween the region of high convective efficiency (the inte- 
rior) and the region of low efficiency (the atmosphere) in 
a much better way than the MLT treatment. 



4. RESULTS AND CONCLUSION 



Results of the calculations of the excitation rates P are 
presented in Fig. (3). The excitation rates P inferred from 
the observations by Chaplin et al (1998) are also repre- 
sented for comparison. 
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Figure 3. Solar p mode excitation rates, P, versus the 
mode frequencies n. Computed P are represented for the 
KMLT model ( dashed line ), the EMLT model (dot dashed 
line ), the KCGM model ( solid line ) and the ECGM model 
(dot dot dashed line). The filled circles with associated 
error bars represent the 'observed' solar values of P{v) 
derived from the a mplitudes and line wid ths of the radial 
modes measured bv \Chaplin et al.\ 111 99^1 



For the MLT treatment, the rates P do significantly de- 
pend on the properties of the atmosphere. Indeed dif- 
ferences in P between the EMLT model and the KMLT 
are found very large. On the other hand, for the CGM 
treatment, differences in P between the ECGM and the 
KCGM models are very small compared to the error bars 
attached to the seismic measurements. 

For the EMLT and KMLT models, P are significantly 
under-estimated compared with the solar seismic con- 
straints obtained from Chaplin et al. (1998)'s measure- 
ments. KCGM and the ECGM models yield values for 
P closer to the seismic data than the EMLT and KMLT 
models. 

We conclude that the solar p-mode excitation rates pro- 
vide valuable constraints on the properties of the super- 
adiabatic region. According to the present investiga- 
tion (focused on local approaches), they clearly favor the 
CGM treatment with respect to the MLT one. 

The remaining discrepancy between computed and ob- 
served P (Fig. 3) is due to the local assumption in the 
convective treatment. Indeed Samadi et al. (2003b) have 
succeeded in reproducing much better the seismic con- 
strai nts by using con straints from a solar 3D simulation 
fsee lSamadi et alJ2b04. present conference). 
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